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Floating Phase in 2D ANNNI Model 



O 

o 

(N 



B 



o 



o 

p 

in 
o 
1^ 
o 



X 



Anjan Kumar Chandra and Subinay Dasgupta 
Department of Physics, University of Calcutta, 
92 Acharya Prafulla Chandra Road, 
Calcutta 700009, India. 

(Dated: February 1, 2008) 

We investigate whether the floating phase (where the correlation length is infinite and the spin-spin 
correlation decays algebraically with distance) exists in the temperature(T) - frustration parameter 
(k) phase diagram of 2D ANNNI model. To identify this phase, we look for the region where (i) 
finite size effect is prominent and (ii) some relevant physical quantity changes somewhat sharply 
and this change becomes sharper as the system size increases. For k < 0.5, the low temperature 
phase is ferromagnetic and we study energy and magnetization. For k > 0.5, the low temperature 
phase is antiphase and we study energy, layer magnetization, length of domain walls running along 
the direction of frustration, number of domain-intercepts that are of length 2 along the direction of 
frustration, and the number of domain walls that do not touch the upper and/or lower boundary. 
In agreement with some previous studies, our final conclusion is that, the floating phase exists, if at 
all, only along a line. 

PACS numbers: 05.70.Jk, 05.10.Ln, 64.60.Fr 



I. INTRODUCTION 

The two-dimensional Axial Next-Nearest Neighbor 
Ising (ANNNI) model (spin = i) is a square lattice Ising 
model with nearest neighbor ferromagnetic interaction 
along both the axial directions and second neighbor anti- 
ferromagnetic interaction along one axial direction. The 
Hamiltonian is, 



n 



(1) 



where the sites {x,y) runs over a square lattice, the spins 
Sx,y are ±1, J is the nearest- neighbor interaction strength 
and K is a parameter of the model. For positive values 
of K the second-neighbor interaction introduces a frustra- 
tion. This is one of the simplest frustrated classical Ising 
model with a tunable frustration and has been studied 
over a long time ([U, 0, Q for review). The most widely 
studied aspect of this model is the phase diagram in the 
T — K phase space, where T stands for temperature. It is 
easy to prove analytically that [l|, [3, at zero tempera- 
ture, the system is in a ferromagnetic state for k < 0.5, 

and in antiphase (-I--I l-H •■■ along x direction 

and all like spins along y direction) for k > 0.5 with a 
"multiphase" state at k = 0.5. (The multiphase state 
comprises of all possible configurations that have no do- 
main of length 1 along the x direction.) 

From Monte Carlo simulations and approximate an- 
alytic calculations, some early studies had proposed a 
phase diagram (Fig. [T]) consisting of a ferromagnetic 
phase (for k < 0.5) and antiphase (for k > 0.5) at low 
temperature along with a paramagnetic phase at high 
temperature (for all k values). The crucial point is, be- 
tween the ordered and the disordered phases for k > 0.5, 
there may be a so-called "floating" phase characterized 
by (i) a spin-spin correlation that decays algebraically 



with distance and (ii) an incommensurate, continuously 
varying modulation. An approximate analytic treatment 
by Villain and Bak \^ predicted that (for k > 0.5) as 
temperature increases from zero, there is a second order 
Pokrovsky-Talapov type commensurate-incommensurate 
phase transition from antiphase to the floating phase fol- 
lowed by a Kosterlitz-Thouless transition from the float- 
ing to the paramagnetic phase. This leads to a phase 
diagram shown schematically in Fig. 1. Several computa- 
tional studies (see fT, '^) also confirmed such a phase dia- 
gram. Later, Shirahata and Nakamura [5] have measured 
the dynamical exponent by studying the non-equilibrium 
relaxation of order parameter at k = 0.6 and 0.8 and 
concluded that the floating phase exists, if at all, over a 
narrow temperature range only. The central problem in 
this study is that the identification of order parameter is 
ambiguous for k > 0.5, in the sense that it is difficult to 
identify a physical observable that relaxes algebraically 
at the critical temperature. One may note that the an- 
tiphase magnetization 



M. 



<2> 



E 



does not satisfy this criterion [5i]. (Here, s"^^^ is the spin 
distribution for perfect aniphase distribution.) Recently, 
a density matrix renormalization group analysis @ has 
also excluded the presence of any incommensurate phase 
over an extended region. 

The 2D ANNNI model is related to the transverse 
ANNNI chain by Suzuki- Trotter transformation 0, H, 
B This quantum Ising model has also been stud- 

ied widely. Numerical and approximate analytic studies 
0, [m [12] had predicted that the floating phase exists 
in the transverse ANNNI chain over a region, as shown 
in Fig. 1. However, recently we have shown 13] that 
for this model, the floating phase exists over a wide re- 
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FIG. 1: Schematic phase diagram of the two-dimensional 
ANNNI model according to previous studies. 

gion extending from k < 0.5 to k > 0.5. There is thus a 
controversy about the existence of algebraically decaying 
phase in 2D ANNNI model. 

In this article, we shall present Monte Carlo simulation 
of 2D ANNNI model with a view to locating the floating 
phase, if any. For k < 0.5, one can easily identify the 
magnetization as the order parameter and for this case 
we have therefore, measured 

1. internal energy 

2. magnetization. 

For K > 0.5, the ordered phase is "antiphase" and it is 
difficult to identify the order parameter unambiguously. 
In this case, we have measured 

1. internal energy, 

2. layer magnetization (magnetization perpendicular 
to the direction of frustration) 

3. length of domain walls running along the direction 
of frustration 

4. number of domain-intercepts that are of length 2 
over a straight line along the direction of frustration 

5. number of dislocations measured as the number of 
domain walls that do not touch the upper and/or 
lower boundary. 

We shall discuss later (Sec. Ill) the significance of these 
quantities in the context of our work. 

From the measurement of a suitable physical quantity 
Q{t) at a time t, the critical point (or, for that matter, 
the critical region) could be identified from the general 
principle that at the critical point the quantity Q{t) — 
Q{oo) is expected to vanish algebraically as a function 
of time t. While this characteristic is handy for the case 
Q{oo) = 0, it is not usable when Q{oo) ^ 0, as very 
large time simulation is required to measure the quantity 
Q{oo) itself. In such cases, we have utilized two essential 



features of critical phenomena : (i) For a finite lattice, 
the correlation time will diverge as 

Hence, the quantity Q{t) will depend strongly on the 
system size only at the critical temperature Tc even at 
finite values of t. (ii) the equilibrium value Q{oo) will 
undergo a sudden change, which is detectable even for 
small size and becomes more and more drastic as the 
system size increases. 

As mentioned above, the results from non-equilibrium 
relaxation study [5i] and the density matrix renormaliza- 
tion group analysis contradicts the previous studies 
BSllBI as regards the extent of the floating phase. This 
paper confirms the conclusion of the former two studies 
by Monte Carlo simulation. One should note that all 
these studies agree at sufficiently low temperatures. We 
measure some observables that play a crucial role in the 
underlying physics and that have not been analyzed till 
now. 

In Sections II and III we shall present the simulation 
studies for k < 0.5 and k > 0.5 respectively. All the simu- 
lation studies were performed with sequential Metropolis 
algorithm using periodic boundary conditions in X and 
Y directions and the results were averaged over 10 to 
50 realizations. In Section IV we shall study the corre- 
spondence between the 2D ANNNI model and the trans- 
verse ANNNI chain and in Section V present conclusions. 
Our final conclusion is that the divergent correlation time 
exists only along a line. The phase diagram obtained 
is presented in Fig. [21 Everywhere in this communica- 
tion temperature is measured in unit of Tc^^ — 0.44069, 
the critical temperature for nearest-neighbor interaction 
(k = 0). This diagram is in qualitative agreement with 
that obtained by Shirahata and Nakamura Q and De- 
rian, Gendiar and Nishino The small difference be- 
tween our results and those obtained by these authors 
seems to be due to the small (~ 1000 x 1000) size of our 
simulation, compared to that of Shirahata and Nakamura 
(- 6399 X 6400). 

II. SIMULATION STUDIES FOR k < 0.5 
A. Energy 

Before we consider the measurements on the 2D 
ANNNI model itself, let us break off for a discussion on 
critical behavior of energy relaxation in general. Internal 
energy E is always unambiguously defined, in contrast 
to the order parameter, which for some phase transitions 
(like our system for k > 0.5) may not be easy to identify 
and measure. However, since E{oo) ^ 0, it is difficult to 
study the time variation of the quantity E(t) — E{oo), as 
mentioned above. The relaxation of energy has also been 
studied elsewhere (l6l.[T7|. 

We shall now consider the case of k (only nearest- 
neighbor interaction), and obtain the exponent, following 
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FIG. 2: The phase diagram of the two-dimensional ANNNI 
model as obtained from the present study, from non-linear re- 
laxation (NER) ^5] and density matrix renormalization group 
(DMRG) 0. The temperature is measured in units of ri"^ — 
0.44069 as mentioned in the text. 
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FIG. 3: Energy relaxation for a 1000 x 1000 lattice for only 
nearest neighbor interaction i.e. k = 0, starting from an or- 
dered (all up) configuration. The straight line is a guide to 
the eye and fits to 0.82f-'' '^*. 



standard scaling arguments Starting from 

the standard diffusion equation for the (non-conserved) 
order parameter ■;/; 

dj^ _ dF_ 
'dt ^ 



one obtains 



dF _ 1 /^V 
9^ ~ ~f \ dt 



where F is the appropriate free energy and F is a param- 
eter of the model. The identity E = d{PF)/d(3, where E 
is the total internal energy then gives, 



dE 1 / B\ fdiPX 



(2) 



The scaling relation for the order parameter may be writ- 
ten as [T9| . 

^{t, e, L, ^o) = h-^'/''^{h-H, fei/'-e, h-^L) (3) 

where e — (1 — T) (T being the temperature measured in 
units of Ti^"*), L is the linear dimension of the system, h is 
the scaling factor, and i/jq is the initial value of the order 
parameter. Also, j3 and v are the static critical expo- 
nents, and z is the dynamic critical exponent. Choosing 
h — i^/^ and suppressing the unimportant arguments L ( 
—>■ oo) and ipQ, we obtain. 



i;{t, e) - t-l'/''' U(l, 0) + t^^^^e^'il, 0) 



(4) 



where ip' is the derivative of with respect to e. Substi- 
tuting this form in Eq. (2), we have, for large time. 



AE - r 



(5) 



where AE is the energy difference E{t) — E{oo) and a = 
1 — (1 — 2(3) /vz = 0.66. Our simulations confirm this 
value of sigma (Fig. [3|) . 

In the case of ANNNI model {k non-zero but < 0.5) we 
study the energy relaxation starting from ferromagnetic 
state. Our results for k = 0.45 is as follows. The E{t) vs 
t curve for an L x L lattice is the same for L = 700, 1000 
and 1200 at a temperature T =0.40 (in units of tP). 
This happens also for T =0.41 but not for T =0.404, 
where the relaxation behavior depends on L to the largest 
extent (Fig. [J). At L = 700, we could evaluate E{oo) and 
found that E{t)—E{oo) vs t plot shows an algebraic decay 
around T=0.40 but the algebraic region is most extended 
for T = 0.404 (Fig. [5]). It is hence concluded that at k = 
0.45, the critical point is Tc = 0.404 ± 0.002. We could 
not however evaluate E{(X)) for L = 1000 or 1200 as they 
involve too much computational time. Moreover, even at 
T = 0.404, there is some anomaly in the sense that the 
decay of E{t) at L = 1200 is faster (rather than slower) 
than that for L = 700 and 1000 (Fig. 15). A simulation 
over a longer time scale might resolve the anomaly. 

It is interesting to note that at k = 0.45, T = 0.404, 
the exponent a has a value 0.15 ± 0.05. Although it is a 
computationally intensive job to determine a accurately, 
we attempted to study the apparent variation of a with 
K and found that it does not change much till k = 0.4 
and afterwards decrease markedly. 



B. Magnetization 

For K < 0.5, the order parameter is magnetisation 
M, whose equilibrium value is zero at the critical point. 
Hence, an easy method of locating the critical point is to 
investigate where the magnetisation relaxes algebraically. 
It is well-known (Eq.(4)) that at Tc the magnetisation de- 
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FIG. 4: Energy relaxation for n — 0.45 at temperature T 
for L X L square lattice. The numbers at the right margin 
indicate T and L values. It is important to note that for 
T =0.40, 0.405, 0.41 and 0.42, the curves for L =1000 and 
1200 coincide with that of L = 700. 
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FIG. 5: Energy relaxation for k = 0.45 for 700 x 700 lattice. 
The numbers indicate temperature. Note that the linear (al- 
gebraic) region is most prominent for T — 0.404. The alge- 
braic region fits to 0.174""'^^. We could not furnish the curve 
for T = 0.405 as the system takes too long time to equilibrate. 



cays as 

M - t-"' (6) 

where a' = P/vz. For k = 0, the value of a' is 0.05734. 

At K = 0.45 the magnetization is found to relax alge- 
braically only around T — 0.405 ± 0.002, which therefore 
is the critical temperature (Fig. [S]). That the critical 
temperature at k = 0.45 lies between 0.40 and 0.41 is 
also verified by the fact that there is a sudden change 
in the equilibrium value of magnetization as tempera- 
ture increases from 0.40 to 0.41, and that this change be- 
comes more and more sudden as the lattice size increases 
(Fig.©. 

The exponent for magnetization decay turns out to be 
a' = 0.02 ±0.005 at k = 0.45, T = 0.404. As for the case 
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FIG. 6: Relaxation of magnetization aX n = 0.45 at for 1000 x 
1000 lattice. The linear (algebraic) region is most prominent 
for T = 0.404 and fits to 1.08f-° °^^ 
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FIG. 7: Equilibrium value of magnetization at k = 0.45 as a 
function of temperature. 



of energy relaxation, it is a computationally intensive job 
to determine a' accurately. Approximate measurements 
indicate that just like cr, the exponent a' also remains 
more or less the same up to k = 0.4 and starts decreas- 
ing markedly at higher n. Further investigations on the 
apparent variation of a and a' is in progress. 



III. SIMULATION STUDIES FOR k > 0.5 

In this case, it is difficult to identify the order pa- 
rameter, as mentioned above. It can be easily proved 
[ESQ that, at low temperature the system is in a per- 
fectly ordered state with like spins along the Y axis and 
+H pattern repeated along the X direction. The do- 
main walls thus run exactly parallel to the Y axis. Mov- 
ing along the X (Y) axis, one finds domains of length 2 
(i), for an L X L system. As one increases the tempera- 
ture, after some temperature Tc (called the lower critical 
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temperature), domains of length larger than 2 appear 
along the X direction. The domain walls now run not 
always parallel to the Y axis. They start from the lower 
boundary and terminate at the upper, but they often 
take small steps parallel to the X axis. Along the Y axis 
the domains are now sometimes less than L in length. 
Villain and Bak [3| pointed out that the number of do- 
mains that do not touch both boundaries is crucial and 
represents some sort of "dislocation" . This number is al- 
most zero immediately above Tc but suddenly increases 
at some temperature T2. While Tc marks a second or- 
der commensurate-incommensurate (Pokrovsky-Talapov 
type) transition, T2 marks a Kosterlitz-Thouless type 
transition. For T > T2 domains of small size appear in 
X and Y direction and the system is in a paramagnetic 
state. For Tc < T < T2, Villain and Bak claimed that, 
the wavelength of modulation changes continuously with 
temperature and the correlation length is infinity, leading 
to a spin-spin correlation that decays algebraically with 
distance. 

To study the two transitions, one at T — Tc and the 
other at r = T2, we measure several quantities always 
starting the simulation with antiphase as the initial con- 
figuration. 

(1) Internal energy is studied in the same way as done 
for K < 0.5. 

(2) Layer magnetization is defined (following [5]) as the 
magnetization of a chain along the Y axis averaged over 
all such chains : 

1 ^ 

m/ = - ^ I m^, I (7) 

x=l 

where, 

1 ^ 

Obviously, this quantity should be 1 for T < Tc and zero 
above T2. Hence, the size dependence of the relaxation 
of mi and the algebraic nature of relaxation of mi would 
indicate a diverging correlation length. Shirahata and 
Nakamura Q have studied this quantity to identify the 
upper transition temperature T2 at k = 0.60 and 0.80. 

(3) The sum of length of the segments of domain walls 
that are parallel to the X axis, divided by the system 
size gives a quantity (say, d^) which is strictly zero in the 
perfect antiphase, but increases suddenly to some non- 
zero value at Tc. Its measurement leads to an estimate 
oiTc. 

(4) Moving in the X-direction, one may note the 
lengths of domain intercepts encountered, and compute 
the ratio 

/2 = n2/nt 

where n2 is the number of domains of length 2 and rit is 
the total number of domains. This ratio is 1 for T < Tc 



and decreases for higher T. We note the region over 
which the relaxation of /2 depends on size or is algebraic 
in nature. 

(5) The spin domains are identified and the fraction 
ifd) of domains that do not touch the upper and/or lower 
boundary is counted. (Those which miss any one bound- 
ary is counted with weightage 1 and those which miss 
both the boundaries are given weightage 2.) This frac- 
tion measures the number of "dislocations" that drives 
the Kosterlitz-Thouless phase transition at T2. We do 
not study the relaxation of this quantity for computa- 
tional limitations, but obtain from simulation the equi- 
librium value away from Tc- Such measurement should 
lead to an estimation of the Kosterlitz-Thouless transi- 
tion temperature T2, if any. 

The study of all these quantities leads us to the con- 
clusion that the floating phase, i.e. the region between 
Tc and T2 extends, if at all, over a temperature range 
less than 0.02. Thus, the region of diverging correlation 
length exists only along a line, up to the accuracy of 
this study. Our study was performed at k = 0.55, 0.60, 
1.0. The resulting phase diagram is presented in (Fig. [2]). 
However, it is interesting to observe that while the width 
of the critical region is 0.02 for k > 0.5, it is much less, 
about 0.001, for k < 0.5. 



A. Energy 

The study of energy relaxation for k > 0.5 follows 
closely the procedure for k < 0.5, with the only differ- 
ence that the initial configuration is now the antiphase. 
For K ~ 0.60, the energy relaxation depends on size pre- 
dominantly at r = 0.43 and 0.44 (Fig. [S]) and for L = 
700, the energy difference E{t) — E{(X)) shows an alge- 
braic decay over an extended region of time at T =0.43 
and 0.44 (Fig. El). Therefore, at k = 0.60 we identify Tc 
as 0.44 ± 0.01. The curves for k = 0.55 are qualitatively 
similar to that for k = 0.60, and Tc could be identified as 
0.37 ± 0.01. An alternative interpretation of the results 
could be that the fioating phase would exist, if at all, be- 
tween T =0.43 to 0.44 (0.37 to 0.38) at k = 0.60 (0.55). 
That the energy does not show any critical behavior over 
an extended range of temperature, seems to indicate that 
the fioating phase does not exist over an extended region. 

We mention that we could not perform the study of 
energy relaxation at k = 1, since for this case a reliable 
data needs averaging over too many configurations. 

The exponent for energy relaxation a' (see Eq.(5)) is 
found to be 1.7±0.1 for k — 0.55 and 0.60. In contrast to 
the findings for k < 0.5, we observe no marked variation 
of a' with K. 



B. Layer Magnetization 

The relaxation of layer magnetization (to/) is qualita- 
tively similar for k = 0.55, 0.60 and 1.00. At k = 0.60 
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FIG. 8: Energy relaxation for k = 0.60 at temperature T 
for L X L square lattice. The numbers at the right margin 
indicate T and L values. For T =0.42 and 0.45, the curves 
for L =1000 and 1200 coincide with that of L = 700. 
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FIG. 10: Relaxation of layer magnetization for k = 0.60 at 
temperature T for L x L square lattice. The numbers at the 
right margin indicate T and L values. 
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FIG. 9: Energy relaxation for k = 0.60 for 700 x 700 lattice. 
The numbers indicate temperature. Note that the linear (al- 
gebraic) region is most prominent for T = 0.44 and fits to 
75 X lO^t"^ 



the relaxation shows critical slowing down and finite size 
effect at T = 0.44±0.01 (Fig.[10l), which is therefore the 
value of T2. The equilibrium value of layer magnetiza- 
tion also shows a sharp fall (that becomes sharper as the 
system size increases) at this temperature, at k = 0.60 
(Fig. [TT|) . For K = 0.55 the value of T2 can be estimated 
in a similar manner to be 0.37 ± 0.01. We could not 
observe the curve for equilibrium value of layer magneti- 
zation at K = 1 because for this case one needs too long 
simulation to get the equilibrium value. 



C. Length of domain walls parallel to the direction 
of frustration 

Like layer magnetisation, this quantity (dx) shows crit- 
ical slowing down and finite size effect at T = 0.44 for 
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FIG. 11: Equilibrium value of layer magnetization at k = 0.60 
as a function of temperature. 



K = 0.60 (Fig. [T^ . The equilibrium value of shows 
a sharp rise at this temperature (Fig. I13p . and this rise 
becomes sharper as the system size increases. Hence, the 
study of dx indicates that is 0.44 ± 0.01 for k = 0.60. 
In the same manner, the value of Tc is estimated to be 
0.37 ± 0.01 for K = 0.55 and 0.76 ± 0.01 for k = 1.00. 



D. Fraction of domains that have length 2 

A study of this quantity (/2) leads to Tc = 0.44 ± 0.01 
for K = 0.60, since size-dependent slowing down of the re- 
laxation of /2 is observed at this temperature (Fig. [T4| . 
Moreover, a sharp fall of the equilibrium value of /2 is ob- 
served at this temperature and this fall becomes slightly 
sharper as the system size increases from 100 to 1000 
(Fig. \TE\i. In a similar manner, the critical temperature 
for K = 0.55 and 1.00 is obtained as 0.37 ± 0.01 and 
0.76 ± 0.01 respectively. 
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FIG. 12: Relaxation ol for k = 0.60 at temperature T 
for L X L square lattice. The numbers at the right margin 
indicate T and L values. For T =0.42, the curves for L —500 
and 1000 coincide. 
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FIG. 13: Equilibrium value of a.t k = 0.60 as a function of 
temperature. 
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FIG. 14: Relaxation of /2 for k — 0.60 at temperature T 
for L X L square lattice. The numbers at the right margin 
indicate T and L values. 
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FIG. 15: Equilibrium value of /2 at «; = 0.60 as a function of 
temperature. 
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FIG. 16: Equilibrium value of /^i at k = 0.60 as a function of 
temperature. 



E. Fraction of domains that do not touch the 
boundary (/d) 

We could not study the relaxation of this quantity since 
averaging over too many reahsations is necessary for a 
reasonably smooth curve. Rather, we could measure the 
equiUbrium value at temperatures where the relaxation 
was not prohibitively slow. It is found that for k = 0.60, 
the critical temperature lies between 0.42 and 0.45, since 
a sudden change is observed in the equilibrium value of 
fd in this temperature range and that this change be- 
comes more and more sudden as the lattice size increases 
fFig. llGp . Similar behaviour is also observed for k = 0.55 
between temperature 0.35 and 0.40. This study could not 
be done for k = 1.0 because of computational limitation 
(one has to average over a good number of realisations to 
get reliable data). Inspite of the computational difficul- 
ties for the study of fd, it is clear that the study of this 
quantity excludes the possibility of critical region with 
width larger than 0.02 along temperature axis. 
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FIG. 17: Schematic phase diagram for Transverse ANNNI 
model (with Hamiltonian Tlq of Eq. (8)), after ref. [13] 



IV. MAPPING TO QUANTUM MODEL 

A quantum Ising model in d dimension is related to a 
classical Ising model in d+1 dimension by Suzuki- Trotter 
transformation and this relation is the basic idea 

behind quantum Monte Carlo algorithm. This transfor- 
mation for the ANNNI model has been discussed in detail 
by Arizmendi [8] and can be summarised as follows. The 
quantum Ising Hamiltonian Ti.q for the one-dimensional 
transverse ANNNI model with N sites is given by, 

N N 

= -J'E(44+i - -14+2) - rE^- (8) 

The ground state for this model is equivalent to a classical 
Ising model in two dimension with Hamiltonian 

mn 
x—1 y—1 



strength of nearest-neighbor interaction in the Y direc- 
tion is p times stronger than that in the X direction, thus 
stabilizing the order in the Y direction. This raises both 
the upper critical temperature T2 (at which the order in 
the Y direction breaks down) and the the lower critical 
temperature Tc (at which the order in the X direction 
breaks down). Both the temperatures are raised by al- 
most the same amount (thus the floating phase always 
has negligible width) but the amount of rise depends on 
the value of p. Now consider the quantum Hamiltonian 
Hq. For a given k if the floating phase extends from Fi 
to r2 (F2 > Fi, the difference F2 — Fi being appreciable), 
then the value of p (say, P2) corresponding to r2 will be 
proportionately larger than the value of p (say, pi ) corre- 
sponding to Fi (see Eq. (11)). Although for every value 
of p, one has Tc ~ T2, the common value of Tc and T2 
at p2 is appreciably larger than that at pi. This shows 
that the presence of floating phase over a wide parame- 
ter region for the quantum model is consistent with the 
absence of the same for the classical model. 



TABLE I: Conclusions regarding the critical region from the 
study of various observables (Sec. II and III here). To is the 
temperature (in units of tJ'^' — 0.44069) at which there is a 
critical region of width < 0.02 





To 


0.45 


0.404 


0.55 


0.37 


0.60 


0.44 


1.0 


0.76 



in the limit (m 00) and (n —^ 00) at a temperature Tq 
where 

ksTq ^ mF 

Jq J' ^ ' 

{ks is the Boltzmann constant) and 

P= ^log[coth(l/m)] (11) 

Obviously, the nearest neighbour interaction in the y di- 
rection is p times the same in the x direction and as 
TO — > 00 this ratio p also tends to infinity. 

As mentioned earlier, the Hamiltonian Hq of Eq. (8), 
describing a transverse ANNNI chain, has also been stud- 
ied widely and numerical and (approximate) analytic 
studies ^ad predicted the existence of float- 

ing phase over a wide region extending from k < 0.5 to 
K > 0.5 (Fig. [T7|) . Why then this extensive presence 
of floating phase does not occur in the phase diagram 
of Fig. [2]? It seems that the explanation is the follow- 
ing. For the classical Hamiltonian Hci of Eq. (9), the 



V. CONCLUSIONS 

(1) For the 2D ANNNI model, diverging correlation 
length and algebraically decaying spin-spin correlation 
exists only along a "line". We present in Table 1 the 
temperature Tq such that within the temperature range 
(To - 0.01) <T <{To + 0.01) lie both the upper critical 
temperature T2 and the lower critical temperature Tc. 
We did not conduct the study at temperature intervals 
smaller than 0.01 except at k = 0.45. 
(2) The phase diagram is topologically different for 2D 
ANNNI model and its Suzuki- Trotter counterpart, the 
ID transverse ANNNI chain. This is in contradiction 
with the fact that often a quantum model and its corre- 
sponding classical counterpart show topologically similar 
phase diagram and have the same critical indices at the 
transition point. Two examples of such behavior are (i) 
2D Ising model and ID transverse Ising model with near- 
est neighbor interaction (k = 0) [lO] and (ii) XYZ chain 
and eight vertex model [2l|. 
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